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Abstract 

We present a description of the equation of state of strongly interacting matter 
within a quasi-particle model. The model is adjusted to lattice QCD data near the 
deconfinement temperature T c . We compare in detail the excess pressure at non- 
vanishing chemical potential and its expansion coefficients with two-flavor lattice 
QCD calculations and outline prospects of the extrapolation to large baryon density. 
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Due to the recent progress of first principle lattice QCD calculations, the 
equation of state (EoS) of strongly interacting matter is now at our disposal 
in some region of temperature T and chemical potential /i. Either the overlap 
improving multi-parameter re-weighting technique [1] or the Taylor expansion 
or hybrids of them [2,3] deliver the pressure, entropy density, quark density, 
susceptibilities etc. The knowledge of these quantities is of primary impor- 
tance for a hydrodynamical description of relativistic heavy-ion collisions, the 
confinement transition in the early universe and possible quark cores in com- 
pact neutron stars. Knowing the phase boundary [2] and the end point of the 
first-order deconfinement transition [4] in the region of non- vanishing chemical 
potential is particularly interesting for the envisaged CBM project at the fu- 
ture accelerator facility FAIR at Darmstadt [5]. In the planned experiments a 
systematic investigation of phenomena of maximum baryon density reachable 
in heavy-ion collisions will be attempted. 

Apart from lattice QCD calculations as purely numerical technique to obtain 
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the EoS, also analytical approaches have been invented to understand the ba- 
sic features. We mention dimensional reduction, resummed HTL scheme, $ 
functional approach, Polyakov loop model etc. (cf. [6] for a recent survey). A 
controlled chain of approximations from full QCD to analytical expressions 
without adjustable parameters describing the lattice data would be of desire. 
Success has been achieved for T > 2.5T C [7] at vanishing baryon density. In 
contrast, the range T > T c , in particular close to T c , is covered by phenomeno- 
logical models [8,9] with parameters adjusted to lattice QCD data at p — 0. It 
is the aim of the present paper to compare in detail the quasi-particle model [8] 
with the recent lattice QCD data [3] in the region around T c by extending the 
focus on finite baryon density. For the first time, we present a quasi-particle 
description of the expansion coefficients of the excess pressure for the quark- 
gluon fluid. Only in such a way an adequate and direct comparison with the 
lattice QCD results [3,10] is possible. 

One way to decompose the EoS is writing for the pressure [2,3] 

p(T,p)=p(T,p = 0) + Ap(T,p), ^^ = f>(^)\ (1) 

p(T,p = 0) was subject of previous lattice QCD calculations (cf. [11] for the 
two-flavor case), while Ap(T, p) became accessible only recently [3,4]. Ap(T, p) 
is easier calculable, therefore, lattice QCD calculations focus on this quantity, 
instead of focussing on p(T,p). In contrast, our model covers p(T,p = 0) and 
Ap(T, p) on equal footing, therefore, we have p(T, p) at our disposal. 

The quasi-particle model of light quarks (q) and gluons (g) is based on the 
expression for the pressure 

P=EPa-B(T^), Pa = ^Jdk^((f+(k) + K(k)) (2) 

where B(T,p) ensures thermodynamic self-consistency [8], s = dp/dT, n q = 
dp/ dp, together with the stationarity condition 8p/5m 2 a = [12]. The k inte- 
grals here and below run from to oo. Explicitly, the entropy density reads 

S = J2a=q,g $a with 1 

1 In massless if 4 theory such a structure of the entropy density emerges by resum- 
ming the super-daisy diagrams in tadpole topology [13], and [14] argues that such 
an ansatz is also valid for QCD. [7] point to more complex structures, but we find 
(2, 3, 4) flexible enough to accommodate the lattice data. Finite width effects are 
studied in [15]. In the $ functional approach the following chain of approximations 
leads to the given ansatz [16]: (i) two-loop approximation for the <E> functional; (ii) 
neglect longitudinal gluon modes and the plasmino branch, both being exponen- 
tially damped; (iii) restore gauge invariance and ultra-violet finiteness by arming 
the self-energies with HTL resummed expressions; (iv) neglect imaginary parts in 
self-energies and Landau damping and approximate suitably the self-energies in the 
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and the net quark number density is 



dkk\f+{k)-f-{k)) (4) 



with degeneracies d q = 12 and d g = 8 as for free partons and distribution 
functions /<f (&) = (exp([E a (k) =F lA/T) + with 5 = +1 (—1) for quarks 
(gluons). The chemical potential is \i for light quarks, while for gluons it is 
zero. 

The quasi-particle dispersion relation is approximated by the asymptotic mass 
shell expression near the light cone 

E 2 a (k) = k 2 + m 2 a , m 2 a (T,fi)=U a (k;T,fi) + (x a T) 2 . (5) 

The essential part is the self-energy n a ; the last term accounts for the masses 
used in the lattice calculation [3], i.e., x q = 0.4 and x g = 0. First direct mea- 
surements of the dispersion relation have been reported in [18]. It should be 
noticed, however, that for the EoS the excitations at momenta k ~ T matter, 
for which more accurate measurements are needed. As suitable parametriza- 
tion of n a , we employ here the HTL self-energies with given explicit T and n 
dependencies as in [8] . The crucial point is to replace the running coupling in 
II a by an effective coupling, G 2 (T,fi). 2 In doing so, non-perturbative effects 
are thought to be accommodated in this effective coupling. This assumption 
needs detailed tests which are presented below. Note that Eqs. (2 - 4) them- 
selves are highly non-perturbative expressions. Expanding them in powers of 
the coupling strength one recovers the first perturbative terms. 

The first expansion coefficients in Eq. (1) follow from (2) as q = "^r~f^rU=o : 



Ci = zJ-J d kk 2 7 - :rv — 2 , (6) 



c 4 
c 6 



3N f 

^ J """" ( e - + 1)5 

4tt 2 T 3 J (e w + l) 4 

f-r / dkk 2 - \e 4uJ - 26e 3uJ + 66e 2uJ - 26e w + 1 (8) 

385tt 2 T 3 J (e w + l) 6 I K ' 



f e 4uJ - 10e 3w + 10e w - l) + 1 f e iuJ - 2e 3w - 6e 2w - 2e w + l) 

3 w v 1 3 uj 2 v 



thermodynamically relevant region k ~ T,fi. The pressure follows by an integration. 
2 As shown in [17], it is the introduced G 2 (T,fi) which allows to describe lattice 
QCD data near T c , while the use of the pure 1-loop or 2-loop perturbative coupling 
together with a more complete description of the plasmon term and Landau damping 
restricts the approach to T > 2.5T C . 
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where a; 
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/T, A 2 = (G 2 /vr 2 + iT 2 d 2 GW)U=o, 



A 4 = (^9 2 G 2 /9/i 2 + ^d 4 G 2 /dp 4 )\ fl=0 . (We have not displayed the terms 
oc x q stemming from the lattice masses; in the calculations presented below, 
however, these terms are included to make the model as analog as possible to 
the lattice performance.) Cj with odd j vanish. In deriving these equations we 
have used the flow equation [8] 



dG 2 
dp 



(It- 
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dT 
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where the lengthy coefficients a^T^riT, p) [16] obey ar(T, p = 0) = and 
a^TiTjfi = 0) = 0. This flow equation follows from a thermodynamic consis- 
tency condition. The meaning of Eq. (9) is to map G 2 , given on some curve 
T(p), e.g., on T(/i = 0), into the p plane to get G 2 (T,p) which is needed to 
calculate p, s, n from Eqs. (2 - 4) at non-vanishing values of p. The terms 
needed in Eqs. (7, 8) follow from the flow equation and its derivatives yielding 
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We adjust G 2 (T) through Eq. (6) to c 2 (T) from [3] for N f = 2. We find as 
convenient parametrization 



G 2 {T) 



GL loop (T), T>T C , 
[G 2 _ loop (T c ) + b(l-T/T c ), T<T C , 

where G2_ loop is the relevant part of the 2-loop coupling 

2ft log(log£ 2 ) 



(12) 
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with fa = (UN C - 2N f )/3, ft = (34iV 2 - 13N f N c + 3N f /N c )/6, and the 
argument C, = X(T—T S ) /T c . T s acts as regulator at T c , and A sets the scale. The 
parameters for N c = 3 are A = 12, T s = 0.87T C , and b = 426.1. Fig. 1 exhibits 
the comparison of Ap and n calculated via Eqs. (2, 4) (dashed curves) or by 
using the expansion coefficients Eqs. (7, 8) (solid curves) with the lattice QCD 



4 




2 



0.8 



1.2 



1.4 

T/T c 



1.6 



1.8 



2 







0.8 



1.2 



1.4 

T/T c 



1.6 



1.8 



2 



Fig. 1. Comparison of the quasi-particle model with lattice QCD results [3] for the 
excess pressure (left panel, for constant n/T) and net quark number density (right 
panel, for constant n/T c ). As for the lattice QCD data (symbols) the quasi-particle 
model results (solid curves) are based on the the expansion coefficients 02,4. For 
comparison, the full quasi-particle model results (dashed curves) are exhibited. 

data [3] based on the coefficients c 2) 4 (symbols). One observes an astonishingly 
good description of the data, even slightly below T c , where the resonance gas 
model [19] is appropriate. 3 Interesting is the deviation of the full model from 
the results based on the truncated expansion in a small interval around T c . It 
should be noted that conceptionally different models [20] reproduce the lattice 
data for Ap and n fairly well above T c . Since for small values of [i the higher 
order coefficients C4 and in particular c§ are less important for Ap and n, a 
more stringent test of the model is accomplished by a direct comparison of 
the individual expansion coefficients q with the corresponding lattice QCD 
results. 

Straightforward evaluation of Eqs. (6 - 8) delivers the results exhibited in 
Fig. 2. Since G 2 (T) was already adjusted to c 2 (T) the agreement is good. It 
should be emphasized that all coefficients q(T) are determined by G 2 {T). That 
means the same G 2 {T) describes also the features of C4 and cq. Particularly 
interesting are the peak of c 4 (left panel of Fig. 2) and the double-peak of 
c 6 /c4 (right panel of Fig. 2) or c 6 (not exhibited) at T c . Numerically, these 
pronounced structures stem from the change of the curvature behavior of 
G 2 (T) at T c which determines the terms d 2 G 2 /<9/i 2 | M=0 and d A G 2 /d/j 4 ]^ via 
Eqs. (9, 10, 11). Neglecting these terms would completely alter the shape of 
c^. Similar to [3], we interpret the peak in C4 as indicator of some critical 
behavior, while the pressure itself is smoothly but rapidly varying at T c . 

In summary we present a quasi-particle model which describes the recent lat- 
tice QCD data for non-vanishing chemical potential remarkably well. Besides 



3 Some reasoning why the model may be applicable also slightly below T c can be 
found in [21]. 
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Fig. 2. The expansion coefficients 02,4 (left panel, data from [3]) and the ratio cq/c^ 
(right panel, data from [10]) as a function of the temperature. 

the excess pressure Ap(T, p) and density n above and even slightly below T c 
at small values of the chemical potential, the individual expansion coefficients 
agree well with the data and turn out to depend on each other. In particular, 
also p(T,p = 0) follows, once G 2 {T) is adjusted. We find a small deviation 
from the data [11] which might be attributed to differences in calculating 
p(T, p — 0) and the coefficients of Ap(T, p) on the lattice. 

Having tested these details of the quasi-particle model, we can directly apply 
the found parametrization and calculate the total pressure at arbitrary baryon 
densities, while lattice QCD calculations are yet constraint to small baryon 
densities. Such applications are of interest for the CBM project at FAIR and 
for studying hot proto-neutron stars and cool neutron stars with quark cores 
and will be reported elsewhere. Another application to cosmic confinement 
dynamics is reported in [22]. These applications need a controlled chiral ex- 
trapolation which must base on improved lattice QCD data. 

Inspiring discussions with F. Karsch, A. Peshier, and K. Redlich are gratefully 
acknowledged. The work is supported by BMBF 06DR121, GSI and EU-I3HP. 
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